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PACS numbers: 75.30.Ds,71.27.+a,75.10.Lp,71.10.Fd 



2 



I. INTRODUCTION 

The role of charge and orbital fluctuations on magnetic couplings and excitations is 
of strong current interest in view of the several zone-boundary anomalies observed in 
spin-wave excitation measurements in the metallic ferromagnetic phase of colossal mag- 
netoresistive manganites." The presence of short-range dynamical orbital fluctuations 
has been suggested in neutron scattering studies of ferromagnetic metallic manganite 
La 1 _ x (Ca 1 _ y Sr y ) x Mn03.- These observations are of crucial importance for a quantitative un- 
derstanding of the carrier-induced spin-spin interactions, magnon excitations, and magnon 
damping, and have highlighted possible limitations of existing theoretical approaches. 

For example, the observed magnon dispersion in the T-X direction shows significant 
softening near the zone boundary, indicating non-Heisenberg behaviour usually modeled by 
including a fourth neighbour interaction term J 4 , and highlighting the limitation of the 
double-exchange model. Similarly, the prediction of magnon-phonon coupling as the origin 
of magnon damping^ and of disorder as the origin of zone-boundary anomalous softening 8 
have been questioned in recent experiments.- - - Furthermore, the dramatic difference in the 
sensitivity of long-wavelength and zone-boundary magnon modes on the density of mobile 
charge carriers has emerged as one of the most puzzling feature. Observed for a finite range 
of carrier concentrations, while the spin stiffness remains almost constant, the anomalous 
softening and broadening of the zone-boundary modes show substantial enhancement with 
increasing hole concentration.-^ 

The role of orbital degeneracy, orbital ordering and orbital correlations on spins dynam- 
ics in ferromagnetic manganites in the entire hole-doping range < x ^ 0.5 is still not 
fully understood due to lack of quantitative comparisons of theoretical calculations with 
experimental results in manganites, as highlighted in the recent review 

Most of the earlier investigations including correlation-induced 0(1/ S) quantum cor- 
rections to the magnon spectrum were carried out in the strong-coupling (double-exchange) 
limit and could not satisfactorily account for the observed anomalous zone-boundary magnon 
softening.- Whereas, with typical parameter values t ~ 0.2—0.5 eV for the mobile e g electrons 
and 2 J ~ 2 eV for their exchange coupling to the localized t 2g core spins, the intermediate- 
coupling regime J ~ W appears to be more appropriate for the ferromagnetic manganties.— 

Recent theoretical studies in the intermediate-coupling regime and including the Coulomb 
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repulsion between band electrons have demonstrated some realistic features such as doping 
dependent asymmetry of the ferromagnetic phase and enhanced magnon softening and non- 
Heisenberg behaviour, thereby highlighting the importance of correlated motion of electrons 
on spin dynamics.— ^2 However, these investigations do not include orbital degeneracy and 
inter-orbital Coulomb interaction, and do not provide any quantitative comparisons with 
experimental results on manganites. Earlier theoretical investigations of the role of orbital- 
lattice fluctuations and correlations on magnetic couplings and excitations have mostly been 
limited to ferro orbital correlations,— and have not addressed the physically relevant stag- 
gered orbital correlations. 

A self-consistent investigation of the interplay between spin and orbital orderings has 
recently carried out using a two-orbital ferromagnetic Kondo lattice model (FKLM) includ- 
ing the Jahn- Teller coupling/ - and doping dependence of the calculated Curie temperature 
was compared with experiments for different ferromagnetic manganites. Finite Jahn- Teller 
distortion and orbital ordering were shown to be self-consistently generated at low doping. 
However, again only ferro orbital ordering was included in this analysis. 

Recent investigations of the correlated motion of electrons using a non-perturbative, 
inverse-degeneracy-expansion based, Goldstone-mode-preserving approach have highlighted 
the role of spin-charge and spin-orbital coupling effects on magnetic couplings and excitations 
in orbitally degenerate metallic ferromagnets.— ~— While spin-charge coupling was indeed 
found to yield strong magnon softening, damping, and non-Heisenberg behaviour, it was 
only on including orbital degeneracy, inter-orbital interaction, and a new class of spin- 
orbital coupling diagrams that low-energy staggered orbital fluctuations, particularly with 
momentum near (tt/2, tt/2, 0) corresponding to CE-type orbital correlations, were found 18 to 
generically yield strong intrinsically non-Heisenberg (1— cos q) 2 magnon self energy correction 
in three dimensions, resulting in no spin stiffness reduction, but strongly suppressed zone- 
boundary magnon energies in the T-X direction. 

In this paper we will apply these new results specifically to the case of ferromagnetic man- 
ganites. We will consider a two-orbital FKLM including intra- and inter-orbital Coulomb 
interactions, with realistic values of hole bandwidth, lattice parameter, and Hund's coupling 
etc. Such a quantitative comparison should provide, within a physically transparent the- 
oretical approach, insight into the role of orbital ordering, orbital degeneracy, and orbital 
correlations on spin dynamics in ferromagnetic manganites, seamlessly covering the entire 
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hole-doping range < x < 0.5 within a single theoretical framework. 

The outline of the paper is as follows. Spin-charge coupling effects on spin dynamics in 
a single-band FKLM are presented in Section II, highlighting the strong differences from 
the canonical double-exchange behaviour. Spin stiffness and spin-wave energies are then 
obtained in the orbitally degenerate state of a two-orbital FKLM in Section III, and com- 
pared with neutron-scattering results for the optimally ferromagnetic manganites (x ~ 0.3). 
Inter-orbital interaction V is included next and Section IV describes spin dynamics in the or- 
bitally ordered ferromagnetic state. The hole-doping (x) dependence of the estimated Curie 
temperature is compared with experimental results for two different cases corresponding to 
wide- and intermediate-band manganites. Finally, the role of (it/2, tc/2, 0)-type orbital cor- 
relations on zone-boundary magnon softening and instability of the FM state near x = 0.5 
are discussed in section V. 



The interplay between itinerant carriers in a partially filled band and localized magnetic 
moments is conventionally studied within the ferromagnetic Kondo lattice model (FKLM): 



involving a local exchange interaction between the localized spins Sj and itinerant electron 
spins <Tj. Due to crystal-field splitting of the degenerate Mn 3d orbitals into e g and t 2g levels, 
the partially filled band in ferromagnetic manganites corresponds to the mobile e g electrons 
whereas the three Hund's-coupled localized t 2g electrons yield the localized magnetic mo- 
ments with spin quantum number S = 3/2. 

We will first consider the single-band FKLM, which is appropriate for the orbitally de- 
generate state of the optimally doped ferromagnetic manganites near hole doping x ~ 0.3. 
As the doped holes are shared equally by the two degenerate orbitals, this hole doping cor- 
responds to band filling n — (1 — x)/2 ~ 0.35 in each band. Orbital ordering and orbital 
correlations will be incorporated later within a two-orbital FKLM including intra-orbital 
and inter-orbital Coulomb interactions. Throughout, we will consider a simple cubic lattice 
with lattice parameter a. 



II. SPIN-CHARGE COUPLING MAGNON SELF ENERGY 
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FIG. 1: The first-order quantum corrections to the irreducible particle-hole propagator </>(q, w). 

The first-order magnon self energy for the single-band FKLM, resulting from quantum 
corrections to the irreducible particle-hole propagator 0(q, u) shown in Fig. 1, was recently 
obtained as:— 
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where refers to the bare magnon energies, e£ = ek — uJS are the exchange-split band 
energies, and the superscripts +/— indicate particle (e^ > cf) and hole (e£ < e?) energies. 
Incorporating correlation-induced self-energy and vertex corrections, the above magnon self 
energy represents a spin-charge coupling between magnons and charge excitations in the 
partially-filled majority-spin band, with the coupling vertex explicitly vanishing for q = 0, 
thus ensuring that the Goldstone mode is explicitly preserved. This result was obtained by 
applying the systematic inverse- degeneracy (1/A/ - ) expansion scheme to a purely fermionic 
representation for the FKLM which allows conventional many-body diagrammatic expan- 
sion. 

The imaginary part of the magnon self energy in the strong-coupling limit: 
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yields finite magnon damping and linewidth at zero temperature, arising from magnon decay 
into intermediate magnon states accompanied with majority-spin charge excitations. On 
the other hand, at the classical level (random phase approximation), magnon damping is 
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FIG. 2: Comparison of the bare (a) and renormalized (b) spin stiffness for the single-band FKLM, 
showing the substantial spin stiffness reduction due to the spin-charge coupling effect. 

possible only due to decay into the Stoner continuum, and is therefore completely absent in 
the intermediate and strong coupling regimes where the magnon spectrum lies well within 
the Stoner gap. 

The small-g behaviour of the magnon self energy in Eq. (2) yields the first-order quantum 
correction to spin-stiffness in d dimensions: 

1 (V6 k ) 2 
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which is down by a factor (1/2S) relative to the bare (classical) spin stiffness: 
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The two competing terms above of order t and t 2 /J correspond to derealization energy 
loss and exchange energy gain upon spin twisting, and determine the overall stability of 
the ferromagnetic state against long wavelength fluctuations. The spin stiffness quantum 
correction is only weakly dependent on J (through uJq), and involves only an exchange 
contribution, which is suppressed for electronic spectral distributions characterized by dom- 
inant saddle-point behaviour (Vek = 0), resulting in enhanced spin stiffness D = —D^ 
and ferromagnetic stability. 

In the double-exchange limit (J — > oo), only the derealization part of the classical spin 
stiffness survives. This contribution has a particle-hole symmetric parabolic band-filling 
dependence. Its behaviour with electron filling n is identical to that with hole doping 
x, vanishing at both ends n = and n = 1 and increasing symmetrically with added 
electrons or holes, leading to the conventional understanding within the DE model that 
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FIG. 3: The bare and renormalized magnon energies and magnon damping for the single-band 
FKLM. The magnon damping linewidths are nearly one-tenth of the renormalized magnon energies 
in the T-X direction. 

ferromagnetism increases with added carriers. However, both the classical and quantum 
exchange contributions involving the (Vek) 2 terms in Eqs. (4) and (5) break this particle- 
hole symmetry, and the stiffness actually vanishes at n < 1 before the band is completely 
filled. 

Since the spin stiffness quantum correction involves only an exchange contribution, in- 
cluding quantum correction is therefore effectively equivalent to enhancing the bare exchange 
contribution by decreasing J. Fig. [2] shows the renormalized spin stiffness D = D^ ' — D^ 1 ' 
in units of ta 2 , evaluated from Eqs. (4) and (5) for different values of J. For band fillings 
below n = 0.5, the renormalized spin stiffness for J ~ W is indeed seen to be close to the 
bare spin stiffness for J/t = 4. 

Figure [3] shows a comparison of the bare and renormalized magnon dispersions in units of 
t, and also the magnon damping obtained for the single-band FKLM. Besides the substantial 
reduction of the renormalized magnon energies, we also find that in the intermediate coupling 
regime (J ~ W), the magnon self energy becomes increasingly non-Heisenberg like, resulting 
in ujx < wr/3, as indeed observed. The magnon linewidth T q obtained from the imaginary 
part of the magnon self energy increases with momentum q, but the ratio r q /u; q is found to 
remain nearly constant at ~ 1/10 upto the zone boundary in the T-X direction, in agreement 
with magnon linewidth measurements.- In the X-M, M-R, and R~r directions, we find the 
ratio r q /w q to be smaller than 0.1. 
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FIG. 4: The renormalized spin stiffness calculated for a two-orbital model, including an AF contri- 
bution due to the Mn super exchange. For J « W, the calculated values are close to the measured 
stiffness ~ 170 meVA 2 for the wide-band compound Lai_ x Sr x Mn03. 

The spin-charge coupling effect on magnon excitations can be readily extended to the two- 
band FKLM involving interaction term — J "Yin Si.(a ia + a^) in Eq. (1). Fig. H] shows the 
renormalized spin stiffness for this two-band model obtained by doubling the spin stiffness 
result D = — from Eqs. (4,5), corresponding to the independent contributions 
of the two degenerate e g orbitals. Here we have taken t = 240meV and lattice parameter 
a = 3.87A for manganites. Corresponding to the AF superexchange interaction between Mn 
spins, an AF contribution to spin stiffness -D^p = zJafScl 2 /6 was also subtracted from the 
above result with Jaf = 2.8meV and lattice coordination z — 6. 

Near optimal band filling n pa 0.35 (hole doping x — 1 — 2n pa 0.3) for ferromagnetic man- 
ganites, the calculated spin stiffness for J pa W is in close agreement with the experimentally 
measured values ~ 170 meVA 2 for the wide-band compound Lai_ x Sr x Mn03. The strong 
suppression near n=0.5 (x=0) seen above with decreasing J could also be a contributing 
factor towards destabilizing ferromagnetism in manganites as x approaches zero, in addition 
to orbital correlations. 

Fig. |5] shows a comparison of the calculated magnon dispersion for the two-band FKLM 
along different symmetry directions in the Brillouin zone with experimental data for the 
narrow-band compound Lao.7Cao.3Mn03 obtained from neutron scattering measurements.- 
The renormalized magnon energy was obtained using w q = 2(u^ — E q ) — zJafS(1 — 7 q ). A 
relatively smaller hopping term t = 180meV (bandwidth W = 12t ~ 2eV) was taken in the 
calculation for this narrow-band compound. Ab-initio calculations also yield an estimated 
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FIG. 5: Comparison of calculated renormalized magnon energy dispersion with experimental data 
for the narrow-band compound Lao.7Cao.3MnO"3 from neutron scattering measurements.— 

bandwidth W = 2eV for La!_ x Ca x Mn03.— The renormalized magnon energies are broadly 
in good agreement with the neutron scattering measurements. Including the zone-boundary 
magnon softening resulting from the spin-orbital coupling effect discussed later will improve 
the agreement near the X point. 
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FIG. 6: Temperature dependence of magnetization calculated using the self-consistent Callen 
scheme with the same set of parameters as above. The calculated Curie temperature is close 
to the measured values for narrow-band manganites. 



From the renormalized magnon energies obtained as above, the finite-temperature spin 
dynamics was investigated using the self-consistent Callen scheme in which the magnetiza- 
tion for a quantum spin-S" ferromagnet:— 

(S - $)(! + $) 25+1 + (S + 1 + $)$ 25+1 
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where the boson (magnon) occupation number: 




TV ^ e^q - 1 



(7) 



q 



in terms of the thermally renormalized magnon energies: 



W q = U}q(S z )/S 



(8) 



Figure |6] shows the temperature dependence of magnetization obtained by self consistently 
solving the coupled set of equations (6)-(8). With the same set of parameters as in Fig. 5 
for the magnon dispersion fit, the calculated Curie temperature (~ 280K) is close to the 
measured value ~ 250K for LCMO. 

As orbital correlations are relatively unimportant at these fillings, the orbitally degenerate 
FKLM calculation presented here provides a good description of ferromagnetic properties. 
Orbital ordering and fluctuations due to inter-orbital interaction and Jahn- Teller distortion 
become important near quarter filling (x ~ 0) and near x ~ 0.5, and therefore must be 
included as both strongly suppress ferromagnetism. Effects of orbital ordering and orbital 
fluctuations on spin dynamics are discussed below. 



Weakly doped manganites such as Lai_ x Sr x Mn03 and Lai_ x Ca x Mn03 exhibit an or- 
bitally ordered state for x 0.2, as inferred from x-ray diffraction and neutron scattering 
experiments.- 1 ^ In this section, we will therefore consider spin dynamics in an orbitally or- 
dered ferromagnet with staggered orbital ordering, and examine the behaviour of the Curie 
temperature with the onset of orbital ordering in the low-doping regime. We therefore 
consider a two-orbital correlated FKLM: 

H = -t^ Aa^w ~ J ^2 S *a*-°v + U n *A*t n *A4 + V n ^ n ^ ( 9 ) 



corresponding to the two e g orbitals fi = a, (3 per site, and include intra-orbital (U) and inter- 
orbital (V) Coulomb interactions. With increasing V, the orbitally degenerate ferromagnetic 
state becomes unstable towards an orbitally ordered ferromagnetic state with staggered 
orbital ordering near half-filling (in the majority-spin band). In the pseudo-spin space of the 
two orbitals, the orbitally-ordered state is exactly analagous to the antiferromagnetic state 



III. ORBITALLY ORDERED FERROMAGNETIC STATE 
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of the Hubbard model with staggered spin ordering.— For simplicity, we consider staggered 
orbital ordering in all three directions, although the end compound LaMnOs exhibits only 
planar staggered orbital ordering. 

The Jahn-Teller-phononic term is also considered to be important in manganites, espe- 
cially in the low and intermediate doping range. 2 - However, the inter-orbital Coulomb inter- 
action has been suggested to be much stronger than the electron-phonon coupling in order to 
account for the observed insulating behaviour in undoped manganites above the Jahn- Teller 
transition and the bond length changes below it.— ^ Generally, Coulombic and Jahn-Teller- 
phononic approaches for manganites have been shown to be qualitatively similar.— Indeed, 
a mean-field treatment of the Jahn- Teller term^ yields an electronic exchange-field term in 
orbital space proportional to the orbital magnetization {ni aa — n^p), exactly as would be 
obtained from the inter-orbital interaction term. 

We consider an orbitally-ordered ferromagnetic state with staggered orbital ordering: 

(n ia ^A = (n m ) B =n + M/2 



where the staggered orbital order M characterizes the density modulation on the two sub- 
lattices A and B. For simplicity, we consider a saturated (half- metallic) ferromagnetic state 
with empty spin-j. bands: (n ia ^) = (n^i) = 0. The ferromagnetic ordering is chosen to be 
in the z direction. 

The effective spin couplings and magnon energies are again determined from the particle- 
hole propagator, now evaluated in the orbitally ordered state. The magnon energies are 
obtained as: 



in terms of the bare particle-hole propagators [x°(q, u>)] for the two orbitals /i, which are 
2x2 matrices in the two-sublattice basis, with [x^aa/bb = \Xo\bb/aa following from the 
orbital-sublattice symmetry in the orbitally-ordered state. The bare particle-hole propaga- 



(n m ) A = (rii a ^)B = n- M/2 



(10) 



u CL = J 2 (25)[A(0)-A(q)] 



(11) 



where A(q) is the maximum eigenvalue of the transverse spin propagator: 




(12) 
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tors are obtained by integrating out the fermions in the orbitally ordered ferromagnetic state 
described by the effective single-particle Hamiltonian matrix: 
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within the self-consistent field (Hartree-Fock) approximation. Here I\ = —JS+Vn and T± = 
JS + (U + V)n are the orbitally independent fields with energy difference T| — T-f- = Un + 2JS 
corresponding to the usual exchange band splitting, and A^ = VM./2 and A^ = (U— V)A4/2 
are the self-consistently determined exchange fields corresponding to the orbital ordering Ai. 
The spin (a) and orbital (fi) indices in the matrix are + and — for spins t and I, and orbitals 
a and (3, respectively. 

The eigenvalues and eigenvectors of the above Hamiltonian matrix yield the bare-level 
band-electron energies and amplitudes for orbital \i and spin a: 



= r. ± v /a I+4 



= - 1 + (JUT- 



= ; 1 



b 2 



(14) 



where © and refer to the two eigenvalue branches (±). Orbital ordering splits the electron 
bands with energy gaps 2A t = VM and 2A^ = (U — V)M for the two spins. At quarter 
filling (n = 1/2), the spin-'f" lower band is completely filled and the upper band © is empty, 
yielding a ferromagnetic insulator, and a metal-insulator transition occurs when orbital 
ordering melts, either driven by band overlap with decreasing V or when the temperature 
exceeds the orbital melting temperature. In our saturated ferromagnetic state with (n^) = 0, 
both branches of the spin- J, electron band are above the Fermi energy. 

The similarity of the orbitally-ordered state with the antiferromagnetic state of the Hub- 
bard model results not only in similar expressions for quasiparticle energies and amplitudes 
as above, but also in an identical self-consistency condition: 

I = v - 

V ^ 
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(15) 



which implicitly determines the magnitude of the orbital order m = 2A^/V as a function of 
the effective interaction strength V. When only nearest-neighbor hopping is present, orbital 
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FIG. 7: Comparison of the doping (x) dependence of calculated T c with experiments for the 
two compounds (a) Lai_ x Ca x Mn03 and (b) Lai_ x Sr x Mn03. Staggered orbital ordering sharply 
suppress ferromagnetism due to band narrowing and reduced electronic derealization. 

ordering at quarter filling sets in for any positive V due to Fermi-surface nesting, whereas a 
finite critical interaction strength is required when frustrating next-nearest-neighbor hopping 
terms are included. 

For finite doping x away from quarter filling, the staggered orbital order parameter m de- 
creases rapidly and vanishes at a critical doping value x c . Although long-range orbital order 
may be susceptible to orbital fluctuations at finite doping in the same way as the staggered 
spin ordering is in the doped antiferromagnet,— this HF result does provide a measure of the 
local orbital order. The pseudogap structure associated with short-range orbital order in the 
metallic phase has been observed in recent high-resolution scanning tunneling microscopy 
and spectroscopy measurements of Lan.7Srn.3Mn03 and Lao.62sCao.375Mn03 thin films.— 

For J ~ W and U ~ W, the Hund's coupling contribution 2JS to the exchange band 
splitting typically dominates over the intra-orbital (Hubbard) interaction contribution Urn. 
Consequently, the calculated magnon energies depend rather weakly on U, and we have 
therefore dropped it for simplicity in the calculations presented below. 

Figure [7] shows a comparison of the doping dependence of the calculated tran- 
sition temperature T c with experiments for the two compounds Lai_ x Ca x Mn03 and 
Lai_ x Sr x Mn03.— ^ Assuming a Heisenberg form o; q = zJS{\ — 7 q ) for the magnon spec- 
trum, the transition temperature T c = (l/3)zJ'S(S + l)/^ q (l — 7q) 1 = (V^Vx was 
calculated approximately from the magnon energy at the X point (it, 0, 0) obtained from 
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Eq. (11). An AF superexchange interaction Jaf between the neighboring Mn core spins was 
also included which correspondingly reduces the magnon energy. Staggered orbital ordering 
was found to persist up to x ~ 0.25, as indicated by the dashed vertical lines, and is clearly 
seen to sharply suppress ferromagnetism due to band narrowing and reduced electronic de- 
localization. In the orbitally degenerate state obtained beyond this critical doping value, 
the calculated Curie temperature is indicated by V = 0. As earlier, we have taken hop- 
ping energies t = 180meV and t = 240meV for the narrow-band and wide-band compounds 
La!_ x Ca x Mn03 and Lai_ x Sr x Mn03, corresponding to bandwidths ~ 2eV and ~ 3eV, re- 
spectively. Other parameters are V/t = 6 and 5, Jaf = 2.8meV and 2.1meV for the two 
compounds. Also, J/t — 4 in both cases, as the bare magnon energy in the orbitally degen- 
erate case then approximately corresponds to the renormalized magnon energy (Fig. 2) for 
J ~ W on including the spin-charge coupling effect as discussed in section II. The anomalous 
zone-boundary magnon softening observed in narrow-band compounds will sharply reduce 
T c as x approaches 0.5. 

The above analogy with the antiferromagnet can be readily extended to inter-orbital 
fluctuations and orbital waves (orbiton) in analogy with transverse-spin fluctuations and 
spin waves.— This analogy then directly yields the orbiton excitation spectrum zJS^l — 7q 
in terms of the orbital exchange coupling J = 4t 2 /V in the strong-coupling limit, quantum 
reduction in the orbital order Ai due to quantum orbital fluctuations, orbital order-disorder 
transition at T c ~ (l/3)zJS(S + 1)/ (1 - 7q)~ 1/2 in the strong-coupling limit (V > t) 
and orbital-order melting induced metal- insulator transition in the weak-coupling limit. 
Similarly, hole motion in an orbitally ordered ferromagnet will lead to scrambling of the 
local orbital order, string of broken orbital bonds, incoherent hole spectral function and 
quasiparticle and band-gap renormalizations due to multiple orbiton emission-absorption 
processes, in analogy with corresponding results for the AF.— Insights into the stability of 
the orbitally-ordered state at finite doping away from quarter filling can similarly be obtained 
from corresponding results for the doped antiferromagnet.— However, if realistic inter-orbital 
hopping terms t a p are included, the continuous orbital-rotation symmetry is broken, resulting 
in a gapped orbiton spectrum, making the orbitally ordered state intrinsically different. 
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IV. CE-TYPE ORBITAL CORRELATIONS AND ANOMALOUS MAGNON 

SOFTENING 

The orbitally degenerate ferromagnetic metallic state of manganites also typically be- 
comes unstable near hole doping x = 0.5 due to onset of CE-type orbital correlations. The 
role of such planar staggered orbital correlations with momentum near (vr/2, tt/2, 0) on spin 
dynamics was recently investigated for a two-orbital Hubbard model: 

H = -t 22{al aa a jaa + al^djpo-) + U 22{n ia fn iai . + n m n in ) + 2j Vn ia n^ (16) 

(ij)er i i 

including intra and inter orbital Coulomb interactions U and V. In the saturated ferromag- 
netic state (band filling n equal to magnetization m), electron correlation effects were incor- 
porated within the framework of a non-perturbative Goldstone-mode-preserving approach.— 
The magnon self energy due to coupling between spin and orbital fluctuations [Fig. 8] was 
shown to generically yield strong intrinsically non-Heisenberg (1 — cosg) 2 momentum depen- 
dence, implying no spin stiffness reduction but strongly suppressed zone-boundary magnon 
energies in the T-X direction. In this section we will extend these results to the two-orbital 
correlated FKLM as in Eq. (9) with similar manganite parameters as used in previous 
sections. 

The correlation-induced spin-orbital coupling magnon self energy for the two-orbital Hub- 
bard model was obtained as: 

Ss P - or Uq)^^ f SP 'to )] l Q -" ^ 

^ ^spin ~r i 'orb 

where f2 or b arid f2 sp i n represent characteristic orbital and spin fluctuation energy scales, the 
angular brackets ( ) refer to averaging over the orbital fluctuation modes Q'. The spin- 
orbital interaction vertex [r sp _ orb (q)], represented diagrammatically in Fig. 8, was shown to 
explicitly vanish at momentum q = in accordance with the Goldstone mode requirement, 
and yields the dominant q dependence of the magnon self energy. 

Although the same overall structure is obtained for the correlated FKLM, there are 
two essential differences. First, the magnon self energy correction for the FKLM is given by 
Smagnon = J 2 (2 S) (j)^ in terms of the correponding quantum correction cf)^ to the irreducible 
particle- hole propagator (instead of U 2 m(j)^ as for the Hubbard model). Second, since the 
bare particle-hole propagator x° f° r the mobile electrons in the FKLM involves the total 
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FIG. 8: The spin-orbital coupling diagrams for the irreducible particle- hole propagator (a) can be 
represented in terms of a spin-orbital interaction vertex r sp _ or t,, the three diagrammatic contri- 
butions to which are shown in (b) involving three-fermion vertices. The missing fourth diagram 
vanishes because of the assumption of complete polarization. 

exchange splitting Um + 2JS, the spin-fluctuation propagator — Ux°) corresponding 

to the U ladders as in Fig. E] yields a gapped spectrum with characteristic gap energy ~ 2 JS 
(instead of the 0(t) magnon energy scale for the Hubbard model). 

Consequently, the first-order spin-orbital coupling magnon self energies for the two models 
are approximately related by: 

CT<.)-^«.)x(^)(^). (18) 

For U = 20t, m « 0.3 (as taken in the Hubbard Model calculation^) and J ~ W, the relative 
factor is about 1/10. For simplicity, we have therefore obtained the spin-orbital coupling 
magnon self energy for the FKLM from the earlier Hubbard model result by multiplying by 
this relative factor. 

Doubling the magnon self energy to account for the independent contributions of the two 
degenerate e g orbitals, the renormalized magnon energy w q = — S sp _ or b(q) for the two- 
orbital correlated FKLM is shown in Fig. [9] for different hole dopings. Here w q is the bare 
magnon energy for the two-band FKLM including the negative contribution z JafS(1 — 7 q ) 
of the superexchange interaction Jaf = 2.1meV between Mn core spins. This bare magnon 
energy was evaluated for J = 4 as it then approximately corresponds to the renormalized 
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FIG. 9: Renormalized magnon energies w q = w q — S sp _ or b(q) including the spin-orbital coupling 
magnon self energy, plotted along the T-X direction for different hole dopings, showing the anoma- 
lous zone-boundary magnon softening due to coupling of spin excitations with CE-type, period 4a, 
planar orbital correlations with momentum modes near (vr/2, tt/2, 0). 

magnon energy (Fig. 2) for J « W on including the spin-charge coupling effect as discussed 
in section II. 

The above analysis clearly shows the importance of CE-type staggered orbital correla- 
tions on the observed anomalous zone-boundary magnon softening. Indeed, a sharp onset 
of CE-type orbital correlations at momentum (vr/2, tt/2, 0) has been observed in neutron 
scattering studies near the Curie temperature,- which gradually diminishes in intensity and 
sharpness with decreasing Ca concentration in Lai_ x (Ca y Sri_ y ) x Mn03 crystals. Interest- 
ingly, this behaviour is exactly similar to the resistivity temperature profile of these crystals 
with varying Ca concentration.— Taken in conjunction with our spin-orbital coupling theo- 
retical result of strong zone-boundary magnon softening produced by such CE-type orbital 
correlations,— this sharp onset of orbital correlations near T c would imply a sudden magnon 
softening and collapse of the ferromagnetic state. As the spin-orbital coupling magnon self 
energy ~ (1 — cosg) 2 does not renormalize the spin stiffness, this behaviour could also ac- 
count for the finite spin stiffness observed just below T c in the low-bandwidth materials, 
which has so far been a puzzling feature. 

In the earlier Hubbard model calculation, the spin-orbital interaction term [r sp _ or b] 2 was 
calculated for the orbital fluctuation mode Q' = (tt/2, tt/2, 0) and an averaging factor of 
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1/10 was included, as obtained on averaging over all orbital fluctuation modes assuming 
a flat distribution. As fluctuation modes near (tt/2, ir/2, 0) become increasingly prominent 
due to onset of CE type combined charge-orbital correlations near half doping (stabilized 
by an inter-site Coulomb interaction Vriirij), the enhanced weight of such staggered orbital 
correlations will substantially increase the magnon self energy and hence the anomalous 
magnon softening as x approaches 0.5. 

V. CONCLUSIONS 

Various aspects of spin dynamics in ferromagnetic manganites were theoretically investi- 
gated in terms of a realistic two-orbital FKLM including intra- and inter-orbital Coulomb 
interactions, staggered orbital ordering, and orbital correlations. For the same set of man- 
ganite parameters, the calculations are in close agreement with recent neutron scattering 
experiments on spin stiffness, magnon dispersion, magnon linewidth, Curie temperature, and 
anomalous magnon softening. The set of manganite parameters taken were lattice parame- 
ter a = 3.87A, hopping terms t = 240meV and t = 180meV (corresponding to bandwidths 
W 3eV and ps 2eV, respectively) for the wide-band (LSMO) and narrow-band (LCMO) 
compounds, Hund's coupling J ps W, inter-orbital Coulomb repulsion V/t ~ 5, and Mn spin 
superexchange interaction Jaf ~ 2.5meV. 

In the orbitally degenerate ground state appropriate for the optimally doped ferromag- 
netic manganites with x ~ 0.3, the spin-charge coupling magnon self energy for the FKLM 
was shown to yield substantial suppression of ferromagnetism in the physically relevant fi- 
nite J regime. Incorporating this magnon energy renormalization, we obtained spin stiffness 
D ~ 150meVA 2 , Curie temperatures T c ~ 300K, and the ratio of the magnon linewidth 
to magnon energy r q /cj q ~ 1/10 for magnon modes upto the zone boundary in the T-X 
direction, and the calculated magnon dispersion fitted well with recent neutron scattering 
data for La!_ x Ca x Mn0 3 . 

The onset of orbital ordering at low doping (x < 0.25) due to the inter-orbital interaction 
V was found to sharply suppress T c due to band narrowing effects associated with the 
staggered orbital order. The doping dependence of the calculated Curie temperature for 
two bandwidth cases corresponding to narrow-band and wide-band compounds was in close 
agreement with the observed behaviour of T c for LCMO and LSMO. 
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Finally, the role of CE-type planar orbital fluctuations with momentum near (ir/2, 7r/2, 0) 
was investigated on spin dynamics in the orbitally degenerate ferromagnet near half doping 
(x ~ 0.5). The magnon self energy due to correlation-induced spin-orbital coupling arising 
from inter-orbital interaction and fluctuation was evaluated for the correlated FKLM from 
the earlier Hubbard model result. For the same set of manganite parameters, the spin-orbital 
coupling magnon self energy with U ~ W was found to yield anomalous zone-boundary 
magnon softening of the same magnitude as observed in neutron scattering experiments. 
Onset of the CE-type charge-orbital correlations and the spin-orbital coupling are therefore 
evidently responsible for the instability of the ferromagnetic state near half doping. 
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